Myriad phases of the Checkerboard Hubbard Model 
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The zero-temperature phase diagram of the checkerboard Hubbard model is obtained in the 
solvable limit in which it consists of weakly coupled square plaquettes. As a function of the on-site 
Coulomb repulsion U and the density of holes per site, x, we demonstrate the existence of at least 
16 distinct phases. For instance, at zero doping, the ground state is a novel d-wave Mott insulator 
(d-Mott), which is not adiabatically continuable to a band insulator; by doping the d-Mott state 
with holes, depending on the magnitude of U , it gives way to a d-wave superconducting state, a 
two-flavor spin-1/2 Fermi liquid (FL), or a spin-3/2 FL. 



The phase diagram of weakly correlated metals tends 
to be relatively simple — superconductivity can occur at 
low temperatures induced by weak attractive interac- 
tions and spin-density waves (SDW) and/or charge den- 
sity waves (CDW) occur under special circumstances 
when the Fermi surface is sufficiently well nested. How- 
ever, with strong interactions, there is no reason not to 
have multiple ordered phases. Inverting this logic, one 
might expect competing phases to be a generic feature 
of strongly interacting systems. So it becomes increas- 
ingly important to find a simple and solvable model of 
strong interactions, which could serve as a paradigmatic 
example showing these complexities of competing orders. 

The Hubbard model is the simplest model of a strongly 
interacting electron gas, but alas, no well-controlled so- 
lution exists in more than one dimension. Here we study 
the Hubbard model, Eq. ([T]), on a checkerboard lattice 
with hopping matrix element t between nearest-neighbor 
sites on elementary square plaquettes and t' between sites 
on neighboring plaquettes. (See Fig. [1]) For t' = t, 
this model reduces to the usual (still unsolved) Hubbard 
model on a square latticcii For t' <^ t, where it is a crys- 
tal of weakly coupled "Hubbard clusters," we are able to 
establish a number of features of the zero temperature 
phase diagram, even for strong interactions, U /t > 1, us- 
ing t' /t as a small parameter. (See Fig. [2]) Particularly 
striking is the large number of zero temperature phases; 
we have established the existence of at least 16 distinct 
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FIG. 1: (Color online) The schematic representation of the 
checkerboard Hubbard model. The hopping amplitudes are 
f = 1 on the solid bonds (blue) and t' <C 1 on the dashed 
bonds (black), r labels sites and R the plaquettes. The lat- 
tice spacing between nearest neighboring sites is set to 1 for 
simplicity. 



phases and there are undoubtedly more in the portions 
of the phase diagram for which we have not yet obtained 
a solution. For instance, at x = and U > 0{t'), the 
ground state is a novel d-wave Mott insulator (d-Mott), 
which is a true new state of matter. More generally, the 
details of the phase diagram depend sensitively on the 
choice of clusters. For example, the phase diagram of the 
dimerized Hubbard modelj^ largely consists of a single 
Fermi liquid (FL) phase, plus a band-insulating (BI) and 
spin-1/2 antiferromagnetic (AF) insulating phase. How- 
ever, multiple competing phases appear to be a common 
feature of the strong interacting limit. 

Recently, there have been a number of cluster- 
dynamical mean-field theory (DMFT) studies of corre- 



1/8 



1/5 



3/10 



3/8 



,FL-1I (1=1/2, 
two bands) 



u 



FL-III (1=3/2) 

IWC + AF 



PS 



WC + AF 



PS 



|af(s=3;2)| 



rf-BCS 



PS 



PS 



WC + AF 



EL-III (s=3/2) 

WC + AF 



FL-II (s=l/2, - 
two bands) 



^AF(s=l/2) I 
^+Qrbilal Ncmatic 



rf-BCS 




rf-BEC 



rf-BEC 



FL 



1/4 

X 



FIG. 2: (Color online) Phase diagram of the checkerboard 
Hubbard model for < a; < 1/2 and aU [/ > 0. Abbre- 
viations: "FL"=Fermi liquid; "s=n/2" = spin-n/2; "PS" — 
phase separation; "AF" = antiferromagnet; "WC" = Wigner 
crystal; "d-BCS" = d-wave superconductor; "d-Mott," "d- 
BEC," and "d-CDW" are phases made of d-wave two-particle 
bound-states (hard-core bosons) which are, respectively, a 
Bose-Mott insulator, a superfluid, and a charge density wave; 
"SDW"= spin-density wave; and "Br'= band insulator. The 
various phases are described in the text. 
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latcd electrons. The complexity of the phase diagram of 
the present model in the t' <^t limit, and the strong de- 
pendence on the precise type of clusters raise questions 
concerning the validity of this approach. Conversely, 
solvable cluster models in the small t' limit can serve 
as interesting benchmark tests for such approximate ap- 
proaches, and for future analog simulations with cold 
fermionic atoms in optical lattices. 

Model Hamiltonian. The checkerboard Hubbard 
model, originally studied in Ref. 0, has Hamiltonian 

^ = - irr'cLcr'. + ^5][rir-l]', (1) 

(rr'>,CT r 

where cj^. creates an electron on site r with spin polar- 
ization (T =t, i and fir = X^cr clcrCra- Here tj-r' is the hop- 
ping matrix element from site r' to r and (rr') denotes 
nearest-neighbor sites, trr' = t ot trr' — t' <^ t when 
(rr') are a pair of sites connected, respectively, by a solid 
bond or a dashed bond shown in Fig. [T] We set t = 1 
as our energy units. Note that the model with uniform 
on-site repulsion U preserves the point group symmetry, 
C4v, of the square lattice. The density of electrons per 
site is defined to be Uei = 1 — x where x is the density of 
"doped holes" per site. Since the model is particle-hole 
symmetric, wc restrict our discussion to < a; < 1- 

Treating t' as a small parameter permits us to solve 
this model using perturbation theory. In the unperturbed 
t' = problem, the 2D lattice consists of decoupled four- 
site square plaquettes. The Hubbard model of a four- 
site plaquette is exactly solvable. The eigenstates of the 
decoupled 2D system are direct products of the eigen- 
states on each plaquette. For most densities, x, the un- 
perturbed ground-state is degenerate, so we use degen- 
erate (or near-degenerate) perturbation theory to derive 
an effective Hamiltonian in the low energy state space. 

Ground states of an isolated plaquette. The eigenstates 
of a single plaquette can be specified by the number of 
doped holes, Qh, the total and z component of the spin, 
and the familiar orbital labels "s" (even under C4 — i.e., 
90° rotation), "px ± ipy" (phase changed by ±7r/2 under 
C4), and "d" (odd under C4). For a single plaquette, 
the Q/i = ground state is unique for any positive U 
and has d-wave symmetry. In the Qh — 1 sector, the 
plaquette ground state for U < Ut ~ 18.6 has spin 1/2 
and px ± ipy orbital symmetry, i.e., it is four- fold degen- 
erate corresponding to spin polarizations s = ±1/2 and 
orbital "chiralities" t = ±1/2. However, for Qh = 1 
and U > Ut, the ground state is spin 3/2 and orbital s 
wave, so it is still four-fold degenerate. When there are 
two holes {Qfi ~ 2), the ground state is unique and has 
s-wave symmetry for all U. The Qh — 3 ground state is 
(trivially) spin 1/2 and s wave. 

In adding holes to the system, the issue arises whether 
it is energetically cheaper to add two holes to one pla- 
quette or one hole to each of two plaquettes. This 
is determined by the sign of the pair binding energy'^ 
A = i;o(0) + £'o(2)-2£'o(l), where i;o((5h) is the ground 



state energy of one plaquette with Qh holes. When 
U < Uc ^ 4.6, A is negative which indicates that doped 
two holes prefer to stay in the same plaquette, effectively 
forming a hard-core boson. When U > Uc, A is pos- 
itive; i.e., two holes repel each other. Note, in Fig. [2l 
that the critical values of U at which level crossings oc- 
cur for the isolated plaquette figure prominently in the 
phase diagram of the perturbed system, as well. 

Effective Hamiltonians. Starting from these states, for 
various ranges of x and U, we can derive the effective low 
energy Hamiltonian in powers of t' . Although this pro- 
cedure reduces the number of dynamical degrees of free- 
dom, it still leaves us with a non-trivial many-body prob- 
lem, which is only solvable in certain cases; the phases 
exhibited in Fig. [2] arc those whose existence we have 
established, but there are compelling reasons to expect 
additional phases to exist in the portions of the x-U plane 
that we have not fully analyzed. We will provide more 
details of the analysis, and a discussion of the regions of 
the phase diagram that have only been partially analyzed 
in a future publication.— 

For < [/ <C 0{t'), the interactions are weak, so the 
zeroth order description is in terms of bands (see below), 
and except at x = [where the Fermi surface (FS) is 
nested] and x = 1/2, where there is a BI, we expect a FL 
description to be valid. 

For < a; < 1/2 and O(Vt') <U <Uc- 0{t'), 2x of 
the plaquettes are occupied by a pair of holes, and (1—22;) 
have no holes. Identifying hole pairs as hard-core bosons, 
the effective Hamiltonian is 

^ blb^,+V^'^ J2 P^P^'^ (2) 

<RR'> (RR'> 

where the bosonic creation operator 6^ creates a hole 
pair on plaquette R, pR, = b^^b^ is the number operator, 
and there is an implicit no-double occupancy (hard-core) 
constraint. Because the zero- hole state has d-wave sym- 
metry and the two-hole state has s-wave symmetry, 6r. 
is a charge 2e field which transforms like a d-wave under 
C4<^ Here t^^^ is the effective hopping of bosons and F'l' 
the repulsion between nearest neighbor bosons, both of 
order t'^ . Their explicit dependences on U are somewhat 
complicated, but we have computed them exactlyj^ 

For < X < 1/2 and Uc - 0{t') < U < Uc + 0{t'), 
both singly charge and doubly charged plaquettes occur, 
so the problem maps onto a rather complicated version 
of the Boson- Fermion model, as discussed in Ref. 0- 

For < a; < 1/4 and Uc + 0{t') <U <Ut, the low en- 
ergy states are a mixture of no-hole and one-hole plaque- 
ttes, where the one-hole states are further distinguished 
by two possible total-spin polarizations s = ±1/2 and 
two orbital chiralities r = ±1/2 — i.e., Px ± ipy Conse- 
quently, the effective Hamiltonian is a two-fiavor version 
of the t-J-V model: 

i/(2) = -f(2) '^R,R'/k./R'.3,-. + ^^'■'^ (3) 

(RR'>,s,T 
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where /r^ ^ t creates a fermion with spin polarization s = 
±1/2 and chiraUty t = ±1/2, and there is a no-double 
occupancy constraint which we have left implicit. Here 
0R,R' is +1 (—1) if the effective bond RR' is along the 
X (y) direction. The hopping parameter t*^^) is order of 
t', while ij(2,2) j-gfgj-g terms of order t'^: 

Hi2.2) ^ J(2) J2 Sh • Sr, + J2 "H«R' (4) 

(RR'> (RR') 

(RR'> 

+ E Sr-Sk,[^Xt^'+J';t^t^'+J:^r^r'], 

(RR'> 

where Sr, tr, and nji are spin, pesudo-spin and density 
operators on plaquette R respectively. Strictly speaking, 
there arc additional "pair-hopping" terms, which wc have 
computed but do not display; for x = 1/4, where ijt^.z) 
is the leading term in the effective Hamiltonian, the pair- 
hopping terms vanish. 

For U > Ut, the one- hole ground state of a single 
plaquette has spin-3/2 so the effective Hamiltonian for 
< X < 1/4 is a t-J-V model for spin-3/2 fcrmions 

^(3)^_i(3) y: /l/r'. (5) 

(RR'>:.s 

(RR'> <RR'> 

where is the plaquette fermion creation operator 
on plaquette R with spin polarization s ~ ±1/2, ±3/2 
(as always there is no double occupancy allowed), and 
Sr and riR are corresponding spin and density opera- 
tors on plaquette R. In this case, the effective hopping 
is order of t'^ while and F^^) are order of t' . 
Consequently, this model always occurs in what, for the 
spin-1/2 model, is considered an imphysical limit J^'^\ 

03)»t(3). 

For U > Uc + 0{t') and 1/4 < x < 1/2, the effective 
Hamiltonian is, again, of the form presented in Eqs. ^ 
and Q (for U < Ut) ov Eq. © (for U >Ut), hut with 
different values of the couplings. Moreover, whereas for 
< X < 1/4 the vacuum state is identified with the a; = 
state with 4 electrons per plaquette, so the mean density 
of fermions per plaquette is Ax; for 1/4 < x < 1/2, the 
vacuum state has 2 electrons per plaquette and the mean 
density of fermions is (2 — 4a:;). 

When X > 1/2, the low energy degrees of freedom are 
always plaquettes fermions, each of which is just an or- 
dinary electron. The effective Hamiltonian is the t-J-V 
model, having the same form as Eq. ([5]), but for spin- 
1/2 fermions with different effective parameters t*^^-* , V^^^ , 
and J*^^'. Here, t'^^^ is order of t', while J^**^ and y^^^ are 
order oft''. So t^^) j(4)^ 

Phase diagram. Much of the structure of the phase 
diagram is obvious from the effective Hamiltonian. Wc 
now sketch some less obvious aspects of the analysis. 



Zero doping. Because the zero-hole ground state of a 
single plaquette is unique and there is a finite gap, at 
X = the unperturbed ground state in the limit t' ^ 
is a direct product state and small t' produces only per- 
turbative corrections which do not change any qualita- 
tive properties of the ground state. This is an insulat- 
ing phase with no broken symmetry. However, despite 
the fact that there arc four electrons per unit cell, this 
state is not adiabatically connected to a BI state, since 
it transforms according to a non-trivial representation of 
the point group: from the c?-wave character of the single- 
plaquette wave function, it follows that the many-body 
wave function changes sign under 90° rotation about a 
plaquette center. One can think of this as a Mott insulat- 
ing state with one d-wave boson per plaquette; hence, we 
call it a "c?-Mott" state. In terms of macroscopic observ- 
able properties, this phase has at least two identifying 
features: (i) The pair-field pair-field correlation function 
('^ot'^rX'^R-i'^R+r'T)! although it falls exponentially with 
distance |R|, has an asymptotic d-wave symmetry (for 
large |R|) upon 90° rotation of r or r' separately, (ii) 
It is an orbital paramagneti^ This d-Mott phase is a 
genuine new state of matter, in contrast with the similar 
state in ladder systems^i^ where there is no C4 symmetry 
to unambiguously distinguish it from a band insulator. 

The fact that the d-Mott phase is not adiabatically 
related to a BI implies that, even for x = 0, there must 
be a phase transition as a function of decreasing U . For 
fixed, small t', when U gets small enough, the gap in 
the isolated plaquette is no longer large compared to t' . 
Specifically, when U <^t' the kinetic energy is dominant 
and the U term can be treated through a weak-coupling 
Hartrec-Fock approximation. Since there arc four sites 
per unit cell, there arc for C/ = four bands as follows: 



ei, = ±y/{t-t')^ +4tt' cos2fc^±y'(t-t')2 +4tt' COS^fcy. 

The top and bottom bands are well separated from the 
two middle bands ek,± « ±2t'(cos^ fc^, — cos^ ky) by a 
gap of approximately 2t. Particle-hole symmetry fixes 
the Fermi energy at for x = 0, so that the FS coincides 
with the lines cos — i cos ky where the two bands 
touch. Consequently, the FS is perfectly nested, and any 
weak positive U induces an SDW ground state ordering 
at (7r/2,7r/2), in which the FS is fully gapped. 

For X = 1/2 there is an insulating ground-state which 
is smoothly connected to the BI state at J7 = 0. For x = 
1, there are no electrons, which is trivially an insulating 
state. 

The hard-core boson model in Eq. ^ has been stud- 
ied^ extensively numerically, and its T = phase dia- 
gram is known. For most x G (0,1/2) (i.e., for boson 
concentration between and 1), it has a uniform super- 
fluid phase, which inherits the d-wave symmetry of the 
bosons, but has no nodal quasiparticles; this is labeled 
d-BEC in Fig. H At x = 1/4, the boson density is 1/2 
per plaquette; in this case, it is easy to see (by map- 
ping the problem to an equivalent XXZ model) that, for 
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y(i) > 2, the ground state is a y^x V2 CDW state of 
d-wave bosons, while for l/'^^^/t'^^^ < 2, the ground state 
is superfluid. At the critical point separating these two 
phases, F*-^-'/^^^'' = 2, the effective Hamiltonian has an 
emergent SU{2) symmetry. In the present case, we find 
that at U = Uc and U ^ Us ~ 2.7, V^^^t'^^^ = 2, and 
that V^(i)/t(i) > 2 for [/,<[/< J7c and y^^V*^^^ < 2 for 
U < Us- Around the d-CDW line in the phase diagram 
shown in Fig.[2l there is a small two-phase coexistence or 
phase separation (PS) region because the transition from 
the d-CDW state to the d-BEC is first-order .liS 

The boson-fermion model which applies for U ^ Uc 
and < a; < 1/2 is quite complicated, and has not been 
fully analyzed. However, as was shown in Ref. 0, it can 
be analyzed for a; ^ taking advantage of the fact that 
the fermions and bosons are dilute. A similar analysis 
applies in terms of the particle-hole transformed model 
for X < 1/2. By increasing [/, there exists a crossover 
(the dotted line in Fig. ^ from the d-BEC phase to a 
BCS-like superconducting phase (d-BCS without nodal 
quasiparticles) . For even larger U , a phase transition into 
a spin-1/2 FL phase with two flavors (bands) of fermions 
(FL-II) is expected for the effective model in Eq. dj) at 
small X. 

The spin-3/2 t-J-V model, Eq. ([5|), is similar to the 
spin-1/2 model in the weak hopping limit; this model 
was studied in detaiL ^^'^^ The ground state depends on 
the ratio of J^^^/F^^-*. In the present problem, we ob- 
tain values of V'^^^ an order of magnitude larger than 
J(3). Thus, the ground state is a spin-3/2 FL (FL-III) 
for X < 1/8. For x ~ 1/% (one fermion on every second 
plaquette), the system forms a -\/2x-\/2 Wigner crystalline 
(WC) state, on top of which the residual antifcrromag- 
netic interactions induce 2x2 AF order (-\/2 x-\/2 WC + 
AF). For X = 1/4, every plaquette is occupied by a single 
fermion, whose spins order to yield a ^/2 x^/2 AF. When 
X — 1/5, there is a concentration 4(1/4 — 1/5) = 1/5 of 
unoccupied plaquettcs which order in a V5 x^/5 WC m 
the background of the ^/2 x ^/2 antifcrromagnetic order 
(•\/5 X y/b WC 4- AF). A particle-hole transformation re- 
sults in a second copy of the same sequence of phases 
(with somewhat different energetics) forl/4<a;<l/2. 



For Uc < U < Ut and x = 1/4, there is one fermion 
per plaquette, so the only terms in the effective Hamilto- 
nian that operate arc those in H^'^''^\ This is a complex 
model with a spin and pseudospin on each plaquette. 
We have solved it by approximating the ground-state as 
a direct product of spin and pseudospin factors. Since 
the spin interactions are AF and isotropic, they form the 
well- understood Ncel ground state of the spin-1/2 AF. in 
which (Sr • Sr') ~ — 1/3 for nearest-neighboring R and 
R'. Then, the effective psueduo-spin Hamiltonian is 

(RR'> a=x,y,z 

where = Ja + JLi^n • Sr-). The ordering of the 
pseudo-spins is determined by the JTq, with largest abso- 
lute value. When U <Un~ 7.3, ~J,^ > -Jy > Jz > 0, 
which indicates that the ground state of pseudospins is 
fcrromagnctically ordered along the x direction. This or- 
bital ordering corresponds to the fact that the electron 
density spontaneously breaks the C4 rotational symme- 
try to C2 with no breaking of translational symmetry, 
so this is an "electron nematic" or " orbital nematic" 
phase.— While for [/„<[/ < Ut, there is no such ne- 
matic ordering.— 

For x = 1/2, with two electrons per plaquette, the 
insulating ground state is adiabatically connected to the 
BI state at [/ = 0. Since the plaquette fermion hopping is 
the dominant term, for 1/2 < x < 3/4 and 3/4 < a; < 1, 
the ground state [except, probably, for a narrow region 
\x - 3/4| < 0{t')] is a spin-1/2 FL, while at a: = 3/4, the 
ground state is a spin-1/2 AF. 

Finite temperature. The finite T phase diagram is also 
interesting^^ and worth future study. For instance, at 
x = 1/4 and for Uc < U < Un, the T = phase is a 
spin-1/2 AF and orbital nematic. At any finite T, the 
spin order is lost, leaving a pure (Ising) nematic phase 
up to a nonzero critical temperature. 
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